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ABSTRACT 

We have reconstructed the spectrum of the afterglow of GRB 970508, on 
May 21.0 UT (12.1 days after the GRB), on the basis of observations spanning 
the X-ray to radio range. The low-frequency power law index of the spectrum, 
a = 0.44 ± 0.07 (F u oc u a ), is in agreement with the expected value a = 1/3 
for optically thin synchrotron radiation. The 1.4 GHz emission is self-absorbed. 
We infer constraints on the break frequecies v c and u m on May 21.0 UT from 
a spectral transition from F v ~ z/~ 6 to F v ~ z/~ L1 in the optical passband 
around 1.4 days. A model of an adiabatically expanding blast wave, emitting 
synchrotron radiation, in which a significant fraction of the electrons cool rapidly 
provides a successful and consistent description of the afterglow observations 
over nine decades in frequency, ranging in time from trigger until several months 
later. 

Subject headings: gamma rays: bursts — gamma rays: individual (GRB 970508) 

1. Introduction 

GRB970508 was a moderately bright 7-ray burst (Costa et al. 1997; Kouveliotou et al. 
1997). It was detected on May 8.904 UT with the Gamma- Ray Burst Monitor (GRBM; 
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Frontera et al. 1991), and with the Wide Field Cameras (WFCs; Jager et al. 1995) on board 
the Italian-Dutch X-ray observatory BeppoSAX (Piro et al. 1995). Optical observations of 
the WFC error box (Heise et al. 1997), made on May 9 and 10, revealed a variable object at 
RA = 06 h 53 m 49^2, Dec = +79°16'19"(J2000), which showed an increase by ~1 mag in the 
V band (Bond 1997), and whose spectrum indicated that the distance to the GRB source 
corresponds to a redshift of at least 0.83 (Metzger et al. 1997). The BeppoSAX Narrow 
Field Instruments (NFIs) revealed an X-ray transient (Piro et al. 1997a) whose position is 
consistent with that of the optical variable, and Frail et al. (1997a,b) found the first GRB 
radio afterglow for GRB 970508; the radio source position coincides with that of the optical 
source (Bond 1997). 

The global properties of the X-ray, optical and radio afterglows of GRBs have been 
successfully described by relativistic blast wave models (e.g. Wijers, Rees and Meszaros 
1997). In this paper we show that the detailed properties of the spectrum of the afterglow 
of GRB 970508, and its evolution, are well described by a blast wave model in which 
significant cooling of the electrons plays an important role (Sari, Piran and Narayan 1998). 
This model predicts the occurrence of an extra break at v c in the spectrum, corresponding 
to the Lorentz factor 7 C above which the electron cooling time is shorter than the expansion 
time of the blast wave. 



2. The X-ray to radio spectrum on May 21.0 UT 

In the second week after the GRB 970508 event many observations of its afterglow 
were made by ground-based and orbiting instruments, which allows us to derive the X-ray 
to radio afterglow spectrum on May 21.0 UT (12.1 days after the event) with only a few 
assumptions. 

During the first month the 4.86 and 8.46 GHz observations by Frail et al. (1997a) 
show strong fluctuations attributed to interstellar scintillation (ISS; see e.g. Goodman 
1997). At 8.46 GHz the interval of 8 to 14 days after the event corresponds to a flare. It 
would therefore appear that the 'true' underlying 4.86 and 8.46 GHz flux densities are 
better represented by long-term average values. To discern, apart from an average flux, a 
possible trend in the 8.46 and 4.86 GHz data of Frail et al. (1997a) we have fitted a power 
law (in time t, measured in days since the burst) to the data for t < 30 days. We find 
F4.86GHZ = (101 ±17) xt(°- 61±0 - 06 ) ^Jy ( Xr 2 cd = 254/10) and F 8 . 46GHz = (321 ± 18) x t (o.274±o.o20) 
yuJy (Xred = 521/25) and we have used these expressions to derive a flux estimate for 
May 21.0 UT at these frequencies. At lower frequencies the afterglow was not (or barely) 
detected during this week. The average Westerbork Synthesis Radio Telesecop (WSRT) 1.4 
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GHz flux density (May 9. 7- Jul 16.44 UT) is < 80 /iJy (2a; Galama et al. 1998a; we note, 
however, that Frail et al. 1997a report a flux density of ~ 100 /iJy). 

Pooley and Green (1997) observed GRB 970508 on six consecutive days (May 17-22) 
with the Ryle telescope at 15 GHz. The source was detected on May 16.52 UT at 1.57 ± 
0.25 mJy, but on none of the other days was an individually significant detection made. 
However, averaging all observations obtained in the interval May 17-22 (8-13 days after the 
event), the source was detected at 660 ± 110 /xJy Again, since ISS can modulate the flux 
densities strongly below~ 20 GHz at this galactic latitude, we assume that the detection 
on May 16.52 UT was boosted by this effect, and that the source-intrinsic 15 GHz flux on 
May 21.0 is better represented by the five-day average stated above. 

Bremer et al. (1998) detected GRB 970508 with the IRAM Plateau de Bure 
Interferometer (PdBI) at 86 GHz but not at 232 GHz. No significant variation is observed 
between the three 86 GHz PdBI detections (May 19.07-May 22.99; center epoch May 21.0 
UT) and, combining them we obtain an average -F 86GHz = 1.71 ± 0.26 mJy. For the 232 
GHz upper limit we use the May 21.15 (day 12.2) upper limit of -F232GHZ < 4.22 mJy (2a; 
Bremer et al. 1998). 

In an observation with the Infrared Space Observatory (ISO) on May 21 the source 
was not detected (Hanlon et al. 1998): preliminary upper limit -F 12Mm < 120/xJy. 

We have fitted a power law, F K = 120^ xr a89±a27 /iJy ( x 2 Icd = 0.23/1), to the three 
K band (2.2 fim) detections of Chary et al. (1998), and find F K = 13^ 5 ^Jy for May 21.0 
UT. 

For an estimate of the optical flux density we used the fit by Galama et al. (1998b) 
to the power-law decay of the differential Coussins R (R c ; A c = 6400 A) light curve, i.e., 
-Fkc — (82.7 ± 2.3) x ^-!- 141 ±o.oi4 ^jy. us j n g an absolute calibration uncertainty of 0.1 mag 
we find 4.8 ± 0.5 /iJy. 

X-ray fluxes (2-10 keV) of the afterglow, obtained with BeppoSAX, have been reported 
by Piro et al. (1997b). These observations were made during the burst, and between 6 
hours and 6 days afterwards. We have inferred the X-ray flux Fx = (7.3 ± 4)xl0~ 4 /iJy for 
May 21.0 UT by extrapolating the power law fit (F x oc t" 1 - 1 * - 1 ) of Piro et al. (1997b). 



3. Discussion 
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3.1. The radio to X-ray spectrum 

The radio to X-ray spectrum on May 21.0 UT is shown in Fig. [I| We have divided 
the spectrum into four parts, corresponding to: the lowest frequencies (I; < 2.5 x 10 9 Hz), 
the low frequencies (II; 2.5 xlO 9 - 10 11 Hz), the intermediate frequencies (III; 10 11 - 10 14 
Hz) and the high frequencies (IV; > 10 14 Hz). Region I has a much steeper spectral slope, 
ai.4-4.MSGH* > 1-1, than region II, a 4 . 8 6~86GH z = 0.44 ± 0.07 (x r 2 cd = 1.4/2). The 1.4 GHz 
flux density (region I) is consistent with the expectation for self-absorbed emission (F u oc z/ 2 ; 
Katz and Piran 1997), i.e. the self absorption break i/ a ~ 2.5 GHz. The low-frequency region 
(II) is in agreement with the expected low- frequency tail of synchroton radiation (F v oc z/ 1//3 ; 
Rybicki and Lightman 1979). The high frequency X-ray to optical slope (a = —1.12 ± 0.07) 
is consistent (Galama et al. 1998a; from here on Paper I) with synchroton radiation from 
electrons, of which a significant fraction cools rapidly, with a power-law distribution of 
Lorentz factors, N(j e ) oc % p (Sari et al. 1998). Optical multi-colour photometry, obtained 
in the first 5 days after the event, showed that the optical spectrum was well represented by 
a power law, and approached a = —1.11 ±0.06 (2.1 to 5.0 days after the GRB; Xred = 3.0/4) 
after reaching maximum optical light (Galama et al. 1998b; we used t > 2.1 days), i.e., 
consistent with the X-ray to optical slope on May 21.0 UT, suggesting that the slope 
remained constant over this period. We note that the extrapolation of the X-ray flux is 
uncertain, but none of the further discussion depends on the validity of the extrapolation, 
since both spectral breaks we infer lie below the optical. In the spectrum the local optical 
spectral slope derived from Galama et al. (1998b) is also indicated. 

For region III, only upper limits are available. Extrapolating from the adjacent regions 
and assuming one single spectral break, we find a peak flux of 4.5 mJy at 10 12 Hz. However, 
the temporal development of GRB 970508 indicates that a second break frequency exists. 



3.2. Evidence for a second spectral break 

As noted before by Katz, Piran and Sari (1998), during the first day the optical to 
X-ray spectral slope of GRB970508 a ~ —0.5; which these authors attribute to rapid 
electron cooling. This is consistent with the UV excess reported by Castro-Tirado et al. 
(1998), and the reddening of the optical spectrum during the first five days (Galama et al. 
1998b). From the data of Galama et al. (1998b) we determine (using Ay < 0.01, based on 
the IRAS 100/iin cirus flux, instead of the value Ay = 0.08 used by Galama et al. 1998b) 
a = —0.54 ± 0.14 for t between and 1.5 days, and a = —1.12 ± 0.04 for t between 1.5 and 
5 days. This suggests that a spectral break moved through the optical passband between 
1.0 and 1.8 days after the burst, separating a range with slope a ~ —0.5 from one with 
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slope a ~ —1.1 (see Fig. |2] upper panel). By May 21 we expect this break to have moved 
downward in frequency to ~ 10 14 Hz (see Fig. 1 and Sect. |3.4| ). Additional observational 
evidence for rapid cooling of a significant fraction of the electrons is given in Paper I. 



3.3. Rapid electron cooling 

Sari et al. (1998) argued that the most energetic electrons may lose a significant 
fraction of their energy to radiation. Their model predicts the occurrence of an extra break 
at v c (corresponding to the critical Lorentz factor 7 C above which cooling by synchroton 
radiation is significant), in addition to the regular break frequency v m (the frequency which 
corresponds to the lowest-energy injected electrons; Lorentz factor 7 m ). The evolution in 
time of the GRB afterglow is determined by the evolution of the two break frequencies: 
u c oc t -1 / 2 and v m oc t -3 / 2 (Sari et al. 1998). We assume adiabatic evolution of the GRB 



remnant; see Sect. |3.6| . Since u m is initially greater than u c (Sari et al. 1998), but decays 
more rapidly, there is a time, to, when the two are equal. For late times t > to we have 
u m < v c and the spectrum varies as F u oc z/ - ^ -1 )/ 2 from z/ m up to v c \ above v c it follows 
F v oc v~ p l 2 and below v m it follows the low-frequency tail, F v oc z/ 1 / 3 (Sari et al. 1998). 



3.4. Constraints on the break frequencies v c and u m for May 21.0 UT 

The 86 GHz PdBI observations do not detect the source before 10 days, followed 
by three detections and subsequently only non detections i.e. at 86 GHz the emission 
peaked between ~ 10 and ~ 14 days (Bremer et al. 1998). This implies that the peak flux 
F v>max ~ 1700/iJy at 86 GHz. The decay after the maximum is quite rapid: for F v oc t s we 
find 5 < -1.1 using the 3cr upper limit on May 28.41 (Bremer et al. 1998). If this maximum 
would correspond to the break frequency u c passing 86 GHz then we would expect the 
decay to be much slower, F u oc t -1 / 4 (Sari et al. 1998). If it reflects the passage of u m the 
subsequent decay would go as F v oc 

t 3(i- P )/4 = r o.9 (g ari et al igg8 )_ The observed decay 
rate is in marginal agreement with the passage of the break frequency u m , but definitely 
excludes that v c was passing 86 GHz. We therefore identify the maximum at 86 GHz with 
the passage of v m at t mi 86GHz ~ 12 days, i.e. at the time of the derived spectrum (May 21.0 
UT). 

There is no evidence for a break at u < 86 GHz (the 86 GHz detection is consistent 
with the low-frequency tail; see Fig. [[]). This implies that v c > z/ m and so for v > v m the 
spectrum is predicted to follow F u oc zy _ ( p_1 )/ 2 = z/~ - 6 (we will use p = 2.2, Paper I). Note 



- 6- 



also that therefore t > t at 12.1 days. The intersection of the high frequency extrapolation 
with F v oc z/~ a6 for v > v m gives an estimate for v c (see Fig. [I]); we find v c = 1.6 x 10 14 
Hz at 12.1 days (May 21.0 UT). Independently, we identify the observed optical spectral 
transition between 1.0 and 1.8 days from a = —0.54 ± 0.14 to a = —1.12 ± 0.04 (Sect. 
372] and Fig. H) with the break frequency v c passing the R c band at t Cj R c . For t < t c Rc we 
expect F v oc z/~°- 6 , while for t > t cRc we expect F v oc v~ v l 2 = u~ lA in the optical passband; 
in excellent agreement with the observed transition. Using t c ^ c = 1.4 days (is c oc t -1 ^ 2 ), we 
find v c = 1.6 xlO 14 Hz at t = 12.1 days (May 21.0 UT). This is in remarkable agreement 
with the value derived above from the intersection of F u oc z/ -06 for v > u m in region III 
with the high-frequency extrapolation. We therefore identify the transition observed in 
the optical passband with the passage of the break frequency u c at £ Cj r c ~ 1.4 days. The 
frequencies z/ a , u m and v c are shown in Fig. ([!]); they determine the division of the spectrum 
into the four regions. We note that the upper limit given for the ISO non detection (Hanlon 
et al. 1998) is a conservative number and further analysis may provide a stronger constraint 
on the shape of the spectrum. 



3.5. The near-infrared range 

At 2.2/im (K band) we expect v c to pass at t c> x ~ 16.5 days (u c oc t~ l l 2 ). Before that 
the K band flux density should decay as F v oc t~ a9 (see above). The decay Fx oc £-°- 89±0 - 27 
(Sect. |3~T| ), fitted to the data of Chary et al. (1998), is in good agreement with this. The 
K to R band spectrum should show a gradual transition between 1.4 days (when v c passes 
Rc) and 16.5 days from F v oc v~ os to F u oc v~ l x . Spectral indices «k-r c , derived from the 
K band data of Chary et al. (1998) and the optical R c light curve of Galama et al. (1998b) 
are in good agreement with this: «k-r c = —0.61 ± 0.12 (t = 4.3), «k-r c = —0.69 ± 0.12 
(t = 7.3), and «k-r c = —0.91 ± 0.27 (t = 11.3; see Fig. |2| lower panel). H band (1.6 /mi) 
HST observations (Pian et al. 1998) detect the OT at H = 20.6 ± 0.3 (F H = 6.2 ± 1.5/xJy, 
t = 24.7) and comparison with Galama et al. (1998b) gives «h-r c — —1-18 ± 0.29 (Fig. ^ 
lower panel). This is consistent with the expectation that v c had by then already passed 
the H band. 



3.6. Adiabatic dynamical evolution of the blast wave and the value of t 

The observed transitions t m! 86GHz ~ 12 days and t C) R c ~ 1.4 days imply that t ~ 
0.006 days (500 s). For t > to the dynamical evolution of the remnant is almost certainly 
adiabatic, while before that time it could possibly be radiative (Sari et al. 1998; but 
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see Meszaros, Rees & Wijers 1998). This means that one would expect nearly all of the 
remnant's observed dynamical evolution to be adiabatic. 

However, the maximum flux at 8.46 GHz is F^ max ~ 700yuJy at t m ,8.46GHz ~ 55 days 
(Frail et al. 1997a); this value of t m ,8.46GHz is as expected from the 86 GHz observations, 
for v m oc t~ 3 / 2 . This is less than the peak of 1700 fiJy at 86 GHz and would argue for 
some radiative losses, since in the perfectly adiabatic case the maximum flux should be 
constant with time, while the peak moves to lower frequencies. Also the rather gradual 
early evolution of the radio light curves of GRB 970508 and the observed transition from 
optically thick to thin emission at 1.4 GHz around t ~ 45 days suggest the possibility of 
radiative evolution (Paper I). However, the absence of a break in the smooth power law 
decay of the optical light curve from 2 to 60 days after the burst (Pedersen et al. 1998; 
Castro-Tirado et al. 1998; Sokolov et al. 1998; Galama et al. 1998b) shows that there is no 
important transition in that period. Also the optical spectra and the optical temporal slope 
indicate adiabatic remnant dynamics with u c > u m during this period (Paper I). Taken 
together, these results indicate that some additional ingredient is needed to fully explain 
the radio and mm behaviour; for example, Waxman, Kulkarni and Frail (1998) have argued 
that the transition from ultrarelativistic to mildly relativistic expansion of the blast wave 
may explain these anomalies. 

4. Conclusion 

We have found that on May 21.0 UT the afterglow spectrum of GRB 970508 contains 
all the characteristic parts of a synchrotron spectrum. We infer that the dynamical evolution 
of the remnant is adiabatic (in agreement with Waxman et al. 1998), but that a significant 
fraction of the electrons have synchrotron cooling times shorter than the remnant's 
expansion time. In accordance with the theory outlined by Sari et al. (1998) this produces 
an additional break in the spectrum. The break frequencies z/ m and z/ c , corresponding to the 
minimum electron energy and the energy above which electrons cool rapidly, respectively, 
are equal at a time to which we find is 500 s. This synchrotron spectrum and adiabatic 
dynamics can explain most of the afterglow behaviour over 9 decades in frequency spanning 
the X rays to radio waves, from the first optical to the late radio data. Some aspects of 
the afterglow are not explained by this model: the flux at v m should be constant, but is 
seen to decrease from 12 days when it lies at mm wavelengths to 55 days when it lies in the 
radio; the self-absorption frequency z/ a is predicted to be time-independent, but is seen to 
decrease with time and the early rise of the radio fluxes is slower than expected (see Paper 
I) and finally the decay after maximum at mm wavelengths is perhaps somewhat faster 
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than expected (see paper I for a discussion on the radio light curves). 

In summary, we find that an adiabatically expanding blast wave emiting synchrotron 
radiation, with a significant fraction of rapidly cooling electrons, describes most of the 
afterglow data on GRB 970508 very well. 
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Fig. 1.— The X-ray to radio spectrum of GRB 970508 on May 21.0 UT (12.1 days after the 
event). The fit to the low-frequency part, a4. 86 _ 8 6GHz = 0.44 ± 0.07, is shown as well as the 
extrapolation from X-ray to optical (solid lines). The local optical slope (2.1-5.0 days after 
the event) is indicated by the thick solid line. Also indicated is the extrapolation F v oc u~ 0M 
(lines). Indicated are the rough estimates of the break frequencies z/ a , z/ m and u c for May 
21.0 UT. 
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Fig. 2. — The spectral slope a as a function of the logarithm of time in days after the burst. 
The upper panel is from Galama et al. (1998b) and shows the spectral slope transition in 
the optical passband (slopes determined from U, B, V, R c and I c band observations; we used 
Ay = 0, see Sect. |3.2| ), where data is collected from Galama et al. (1998b; *), Sokolov et 
al. (1997; x), Djorgovski et al. (1997, •), Metzger et al. (1997b, o) and Castro-Tirado et al 
(1998; □). The transition from F v oc z/~ ' 6 to F u oc z/~ L1 (dashed lines) occurs at roughly 1.4 
days. We expect a similar transition in the spectral slope derived from K-R c (lower panel) to 
occur between 1.4 and 16.5 days (see Sect. |3T5| ). The lower panel shows the spectral indices 
«k-r c (I), derived from the K band data of Chary et al. (1998), and «h-r c (A), derived 
from the H band data of Pian et al. (1998), with respect to the the optical R c light curve of 
Galama et al. (1998b). 



